Control of unstable steady states by time-delayed feedback methods 



in 
o 
o 

(N 
< 



5-1 

43 



03 



C 

O 

o 



> 

m 
oo 
o 
in 
o 
-i— > 

03 



o 
o 



X 



P. HoveE and E. Scholl 

Institut fur Theoretische Physik, TV Berlin, Hardenbergstrafle 36, D-10623 Berlin, Germany 

(Dated: February 2, 2008) 

We show that time-delayed feedback methods, which have successfully been used to control un- 
stable periodic ortbits, provide a tool to stabilize unstable steady states. We present an analytical 
investigation of the feedback scheme using the Lambert function and discuss effects of both a low- 
pass filter included in the control loop and non-zero latency times associated with the generation 
and injection of the feedback signal. 
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I. INTRODUCTION 

Starting with the work of Ott, Grebogi and Yorke 
a variety of methods for controlling unstable and chaotic 
systems have been developed in the last fifteen years and 
applied to various real systems in physics, chemistry, biol- 
ogy, and medicine [2j,|3|,|j]. Pyragas |5| introduced a time- 
delayed feedback scheme that stabilizes unstable periodic 
orbits (UPO) embedded in a chaotic attractor by con- 
structing a control force from the difference of the current 
state to the state one period in the past. This method 
is known as time-delay auto synchronization (TDAS) and 
was improved by Socolar et al. [|| by considering multi- 
ple delays in form of an infinite series (extended TDAS or 
ETDAS) or an average of N past iterates (N time delay 
auto synchronization or NTDAS) or coupling matri- 
ces (generalized ETDAS or GETDAS) @. In parallel to 
the control of UPOs, the stabilization of unstable steady 
states (USS) became a field of increasing interest. 

One of the methods to control an USS introduced by 
Bielawski et al. uses the derivative of the current state 
as source of a control force 0. It can be shown, however, 
that this derivative control is sensitive to high frequency 
oscillations [T^j an d thus not robust in the presence of 
noise. Another control scheme was given by calculating 
the difference of the current state to a low-pass filtered 
version [Tlj . 

Although the effects of time-delayed feedback schemes 
on the stability of periodic orbits are understood to large 
extend [13, US LL-J US HH > much less is known in the case 
of a fixed point. There are some results discussing the 
application of the ETDAS control method 0] and nu- 
merical simulations of Chua's circuit 0, but a detailed 
theoretical investigation is still missing. 

The purpose of this paper is the analytical and numer- 
ical study of the TDAS method, which was originally in- 
vented to control unstable periodic orbits [f| , in applica- 
tion to unstable fixed points, including latency and filter- 
ing effects. The paper is organized as follows. In Sec. |n] 
we will introduce the system's equations and the control 
force. In Sec. IIIII we will investigate the domain of con- 



trol in dependence on time delay and feedback gain and 
present analytical solutions of the characteristic equation 
using the Lambert function. Further, we will consider ef- 
fects of non-zero latency times and additional low-passs 
filtering in Sec. HVI and Sec.[V] respectively. 



II. CONTROL BY TIME-DELAYED FEEDBACK 

We consider a general dynamic system given by a vec- 
tor field f: 

x = f(x) (I) 

with an unstable fixed point x* given by f(x*) = 0. The 
stability of this fixed point is obtained by linearizing the 
vector field around x*. Without loss of generality, let 
us assume x* = 0. In the following we will consider the 
generic case of an unstable focus for which the linearized 
equations in center manifold coordinates x, y can be writ- 
ten as 



x = A x + lj y 
V = -ojx + Xy, 



(2) 



where A and uj are positive real numbers. They may be 
viewed as parameters governing the distance from the 
instability threshold, e.g., a Hopf bifurcation of system 
(I), and the intrinsic eigenfrequency, respectively. For 
notational convenience, Eq. can be rewritten as 



x(i) = A x(t). 



(3) 



The eigenvalues Ao of the matrix A are given by Ao = 
A ± iui, so that for A > and w / the fixed point is 
indeed an unstable focus. A vanishing imaginary part, 
i.e., u> = 0, would correspond - in the case of a UPO - 
to an orbit without torsion for which TDAS fails [l2| . 
We note that the same holds for USSs and therefore we 
restrict our investigation to lo 7^ 0. 

We shall now apply time-delayed feedback control [B| 
in order to stabilize this fixed point: 
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x(t) = Xx(t)+uy(t)-K[x(t)-x(t-T)] 
y(t) = -u>x{t)+Xy{t)-K[y(t)-y{t-r) 



(4) 
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FIG. 1: Control of an unstable focus with A = 0.5 and ui = n 
in the configuration space for different values of the feedback 
gain K. Panels (a), (b), (c), and (d) correspond to K = 

0, 0.2, 0.25, and 0.3, respectively. The time delay r of the 
TDAS control scheme is chosen as 1, corresponding to r = 
To/2 = tt/lu. 

where the feedback gain K and the time delay r are real 
numbers. The goal of the control method is to change 
the sign of the real part of the eigenvalue. 

Since the control force applied to the i-th component 
of the system involves only the same component, this 
control scheme is called diagonal coupling which 
is suitable for an analytical treatment. Note that the 
feedback term vanishes if the USS is stabilized since 
x*(t— r) = x*(t) and y*(t—r) = y*(t) for all £, indicating 
the non-invasiveness of the TDAS method. 

Figure ^ depicts the dynamics of the controlled unsta- 
ble focus (A = 0.5 and u> — it) in the x-y plane for differ- 
ent values of the feedback gain K. Panels (a) through (d) 
correspond to increasing K. The time delay of the TDAS 
control scheme is chosen as r = 1 in all panels. Panel (a) 
displays the case of the absence of control, i.e., K = 0, 
and shows that the system is an unstable focus exhibit- 
ing undamped oscillations on a timescale To = 2tt/uj = 2. 
It can be seen from panel (b) that increasing K reduces 
the instability. The system diverges more slowly to in- 
finity indicated by the tighter spiral. Further increase 
of K stops the unstable behavior completely and pro- 
duces periodic motion, i.e., a center [see panel (c)]. The 
amplitude of the orbit depends on the initial conditions, 
which are chosen as x — 0.01 and y = 0.01. For even 
larger feedback gains, the trajectory becomes an inward 
spiral and thus approaches the fixed point, i.e., the focus. 
Hence the TDAS control scheme is successful. 

An exponential ansatz for x(t) and y(t) in Eq. 

1. e., x(t) ~ exp(At), y(t) ~ exp(At), reveals how the 
control force modifies the eigenvalues of the system. The 
characteristic equation becomes 

[A + K (1 - e- Ar ) - A] 2 + lo 2 = 0. (5) 

so that the complex eigenvalues A are given in the pres- 




FIG. 2: (Color online) (a) Largest real part of the complex 
eigenvalues A vs. r for A = 0.5 and ui = ir for different K. 
Some lower eigenvalues are also displayed for K = 0.3 (green 
online), (b) Time series of the x-component of the unstable 
focus: The solid line (red online) corresponds to x(t), the 
dashed line (green online) to the delayed x component x(t—r) 
with t — 1. The parameters of the unstable focus and the 
control scheme are as in panel (d) of Fig. 



ence of a control force by the implicit equation 

X±iuj = A + K(l-e- Ar ). (6) 

Using the Lambert function W, which is defined as the 
inverse function of g(z) — ze z for complex z [2(], 0, [2^, 
HE 01 j Eq- @ can be solved analytically 

At = W ( K KTe- {x±iu)T+KT ^ + (A ± iuo)r - Kt. (7) 

Panel (a) of Fig.[3shows the dependence of the largest 
real part of the complex eigenvalues A upon the time 
delay r according to Eqs. ©, 10 for A = 0.5 and 
u> = 7T. The solid curve corresponds to a feedback gain of 
K = 0.3, the dashed curve to K — 0.25, and the dotted 
curve to K = 0.2. All curves start at Re(A) = A for 
t = 0, i.e., when no control is applied to system. For 
increasing time delay, the real part Re (A) decreases. It 
can be seen in the case of K = 0.3 that there exist values 
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of the time delay for which Re (A) becomes negative, and 
thus the control is successful. The curve for K = 0.25 
shows the threshold case where Re(A) becomes zero for 
t = 1, but does not change sign. The TDAS control 
scheme generates an infinite number of additional eigen- 
modes. The corresponding eigenvalues are the solutions 
of the transcendental Eq. © . The real parts of the eigen- 
values all originate from — oo for r — 0. Some of these 
lower eigenvalues are displayed for K — 0.3. The differ- 
ent branches of the eigenvalue spectrum originate from 
the multiple-leaf structure of the complex Lambert func- 
tion. The real part of each eigenvalue branch exhibits 
a typical nonmonotonic dependence upon r which leads 
to crossover of different branches resulting in an oscilla- 
tory modulation of the largest real part as a function of 
r. Such behavior of the eigenvalue spectrum appears to 
be quite general, and has been found for various delayed 
feedback coupl ing schemes, including the Floquet spec- 
trum of UPOs jjjll2a| and applications to noise-induced 
motion where the fixed point is stable p^ . 

The notch at t = 1 corresponds to Fig. ^ so that at 
this value of r the solid, dashed, and dotted curves cor- 
respond to panels (d), (c), and (b) of Fig. ^ respectively. 
The notches at larger r become less pronounced leading 
to less effective realization of the TDAS control scheme, 
i.e., a smaller or no r-interval with negative Re(A). 

In the case of an UPO the optimal time delay is equal 
to the period of the orbit to be stabilized. Note that in 
the case of an USS, however, the time delay is not so 
obviously related to a parameter of the system. We will 
see in Sec. IIIII which combinations of the feedback gain 
K and the time delay r lead to successful control. 

Panel (b) of Fig. [3 displays the time evolution of x(t) 
and its time delayed counterpart x(t — r) in the case of a 
combination of K — 0.3 and r = 1 that leads to successul 
control as in panel (d) of Fig. The x component of 
the control force can be calculated from the difference 
of the two curves and subsequent multiplication by K. 
Since x(t) tends to zero in the limit of large t (the system 
reaches the focus located at the origin), the control force 
vanishes if the system is stabilized. Thus the control 
scheme is non-invasive. Note that the current signal (red 
online) and its delayed counterpart (green online) are in 
anti-phase. This observation will become important in 
Sec. El 



In the following discussion, it is helpful to consider the 
real and imaginary part of Eq. © separately in order 
to gain some analytic information about the domain of 
control: 



K [1 



-pr 



Ke~ pT sin(gr) 



A 

LO 



(8) 



III. SHAPE OF THE DOMAIN OF CONTROL 

This section is focused on the construction of the do- 
main of control in the K-t plane. The calculation can be 
done analytically for special points by using, for instance, 
that p = at the threshold of control. Furthermore, we 
will present an expansion around the minimal value of K 
that reveals further details of the shape of the domain of 
control. 

At the threshold of control the sign of the real part p 
of the exponent A changes. Therefore setting p to zero 
in the real and imaginary parts, respectively, of Eq. (JHJ 
yields 



and 



A = K [1 - cos(gr)] 



lo = q + ATsin(gr). 



(9) 



(10) 



Since the cosine is bounded between —1 and 1, the 
following inequality follows from Eq. JHJ 



A 



< K. 



(11) 



A/ 2, for which the 
It corresponds to 



Thus a minimal value of K, K m , 
control starts, can be inferred 
qr = (2n + 1)tt for n = 0, 1, 2, . . .. 

In order to express the values of the time delay r that 
correspond to the minimal K in terms of the parameters 
of the uncontrolled system, it is useful to consider even 
and odd multiples of 7r for qr, i.e., qr = 2mr and qr — 

(2n + l)n for n = 0, 1, 2, In both cases, the imaginary 

part of Eq. (JSJ leads to q = lo. Hence, in the latter case, 
the time delay r for K m i n — A/2 becomes 



-(2n+l). 



(12) 



The last expression can be rewritten using the uncon- 
trolled eigenperiod To 



T = Tr 



2n + l 
~ 2~ ' 



where Tq is defined by 



To 



2tt 



(13) 



(14) 



This discussion has shown that K = A/2 and r = 
To(2n + l)/2 with n = 0, 1, 2, . . . correspond to points of 
successful control in the K-t plane with minimal feed- 
back gain. 

For even multiples, i.e., qr — 2nir for n — 0, 1, 2, . . ., 
no control is possible for finite values of since 



K -X 



K 



cos(gr)| 



qr— 2mr 



with A 



p + iq. 



O 1 



K 



(15) 
(16) 
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which cannot be satisfied for A ^ and finite K. Further- 
more, Eq. Ijl U| l yields that for time delays, which are in- 
teger multiples of the eigenperiod, i.e., r = Tqu = 2im/uj 
with n = 0, 1, 2, . . ., the control scheme fails for any feed- 
back gain. Note that this failure appears to be related to 
the case of torsion free UPOs, where it has been shown 
that u> 7^ is a necessary condition for control 

Another result that can be derived from Eq. 10 is a 
shift of q for increasing K. For this, taking the square of 
the real and imaginary part of Eq. JBJ and using trigono- 
metrical identities yields 



q = lo T y/(2K-\)\. 



(17) 



Inserting Eq. (|17|l into the real part of Eq. J^J leads to 
an explicit expression for the dependence of time delay r 
on the feedback gain K at the threshold of stability, i.e., 
the boundary of the control domain p = 0, 



K — X 



K 

r(K) = 



= cos(gr) 



. K , 



■ T VW - A) A 



(18) 
(19) 



In order to visualize the shape of the domain of control 
we will investigate how small deviations e > from K m i ni 
i.e, K = A/2 + e, influence the corresponding values of 
the time delay r. For this, let rj > be small and r = 
^(2n+l)±?7 a small deviation from r at K min . Inserting 
the expression for K and r into Eq. I|18|l yields after some 
Taylor's expansions 



1 



A 



-e = -1 



W7j=F-(2ri+l)V2A v / e (20) 



±— ^ + ^-(2n+l)VA 



Ve- (21) 



This equation describes the shape of the domain of con- 
trol at the threshold of stabilization, i.e., p = 0, near the 
minimum K value at r = Tb(2n+l)/2 in the K-t control 
plane. Small deviations from r at K m i n are influenced 
by the square root of small deviations from the minimum 
feedback gain. 

Figure|3]displays the largest real part of the eigenvalues 
A in dependence on both the feedback gain K and the 
time delay r for lo = ir and two different values of A, and 
summarizes the results of this section. The values of A 
are calculated using the analytic solution of Eq.@). 
The two-dimensional projections at the bottom of each 
plot extract combinations of K and r with negative p, 
i.e., successful control of the system. In the absence of 
a control force, i.e., K = 0, the real part of A starts at 
A. Increasing the feedback gain decreases Re(A). For 
K = K m i n = V^j the real part of the eigenvalue reaches 
for certain time delays, i.e., r = Tb(2n + l)/2 with 
n = 0,1,2,.. ., and then changes sign. Thus, the system 
is stabilized. For values of the feedback gain slighty above 
the minimum value K m i n , the domain of control shows a 



Re(A) 




(b) 



FIG. 3: (Color online) Domain of control in the K-t plane 
and largest real part of the complex eigenvalues A as a func- 
tion of K and r according to Eq. 0. The two-dimensional 
projection at the bottom shows combinations of r and K, for 
which Re(A) is negative and thus the control successful [Panel 
(a): A = 0.5 and u = it, panel (b): A = 0.1 and lo = 7r]. 



square root shape. It can be seen that for time delays of 
r = TqU the largest real part of the eigenvalues remains 
positive for any feedback gain. For a smaller value of A 
(Fig. Eb), i.e., closer to the instability threshold of the 
fixed point, the domains of control become larger. 

An example of the combination of minimal feedback 
gain K m i n = A/2 and corresponding time delay r = 
To(2n + l)/2, n = 0,1,2,... is shown in panel (c) of 
Fig.Hl where K = A/2 = 0.25 and r = T /2 = tt/lu = 1. 
It describes the control threshold case between stable and 
unstable fixed point. 



IV. LATENCY TIME EFFECTS 

In this section we will consider non-zero latency times, 
which can be associated with the generation and injec- 
tion of the feedback signal [2^. It has been shown ex- 
perimentally |2!| in the case of an UPO that latency can 
have important effects on the controllability of the sys- 
tem and might limit the success of the time-delayed feed- 
back method. A theoretical explanation can be found in 
[30l l3l| . Here we will discuss how latency times change 
the domain of control in the case of an USS. 
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FIG. 4: Largest real part of the eigenvalues A vs. r for 
A = 0.5, lo — tt, and K — 0.3 as given by Eq. $23$. The 
solid, dashed, dotted, and dash-dotted curves correspond to 
a latency time of S = 0, 0.1, 0.2, and 0.3, respectively. 



FIG. 5: Minimal feedback gain K vs. relative latency 5/t for 
A = 0.5 and uj — tv according to Eq. 12411 . The shaded areas 
show the domain of control for suitably chosen r. 



The latency time 8 can be included as an additional 
time delay in the control force of Eq. which then be- 
comes 

*<*-«>--* (51: 3 < 22 > 

leading to a characteristic equation similar to Eq. JfjJ) but 
with an additional exponential factor 



X±iuj = A + Ke- Ad (l-e- AT ) (23) 
or, separating into real and imaginary parts, 



p + K 
q-K 



?r pS cos(qS) - e~ piT+S) COS (q(T + 8)) 



A 



' sin(q8) 



_ p -p{r+6) t 



Mr + 8)) 



where p and q denote the real and imaginary part of A, 
respectively. 

Figure 0] displays the dependence of the largest real 
part of the complex eigenvalues A on the time delay r 
according to Eq. ill'MI) for A = 0.5 and u> — it, and different 
values of the latency time 8. The values of the eigenvalues 
are calculated by solving Eq. H23fl numerically. The solid, 
dashed, dotted, and dash-dotted curves correspond to 
(5 = 0, 0.1, 0.2, and 0.3, respectively. The case with zero 
latency time is also displayed; it corresponds to the solid 
curve in Fig. HJa). 

It can be seen that increasing latency time shifts the 
minimum of p to smaller values of r and reduces the 
intervals of r, for which p is negative, i.e., for which the 
control is successful. The dash-dotted curve (8 = 0.3) 
shows a case where no control is possible since the largest 
real part of the complex eigenvalues remains positive for 
all time delays r. 

Solving Eq. I|23|l analytically by using the Lambert 
function, as for Eq. JfJJl, is not possible in the case of 



non-zero latency times due to the additional exponen- 
tial term exp(— AS). In order to understand the effects 
of non-zero latency times on the value of the minimal 
feedback gain, we evaluate the real part of Eq. (|23ll at 
the threshold of control, i.e., p = 0. It was shown in 
Sec. IIIII that K becomes minimal if qr = ir(2n + 1) for 
n = 0, 1,2, . . . [see Eqs. l|§fl.l|TT)l]. This value of qr yields 
for non-zero latency times 



^■min (^) 



A 



2cos[?r(2n+ 1)£] 



> 



A 



(24) 



This shows that non-zero latency times shift the mini- 
mal feedback gain K min , for which the control scheme is 
successful, to larger values. 

The effects of the cosine function in Eq. (|24J) can be 
understood by considering Fig. [2tb) and Fig. |3] which 
shows the dependence of the minimal feedback gain K m i n 
on the latency time for A = 0.5. Increasing latency time 
increases the value of K min . If S becomes larger than 
half the time delay r, K m i n changes sign and control is 
possible only for negative K with K < K m i n and suitably 
chosen r. Note that in Fig. EJtlie difference x(t) — x(t — r) 
has to be positive for successful control. For 0.5r < 6 < 
1.5r both x(t — 6) and x(t — S — r) becomes closer to zero. 
Therefore in order to achieve control, the feedback gain 
becomes larger. In the limit 8/t — ► 1/2 the difference 
x(t — 5) — x(t — t — 5) vanishes and thus the minimal 
feedback gain diverges. For even larger values of 8 the 
above-mentioned difference changes its sign forcing K m i n 
to do the same. Otherwise the control scheme would 
generate a force that pulls the system away from the 
target fixed point. 

Figure[fj]depicts the domain of control for latency times 
of 8 = 0.1 [panels (a) and (c)] and 8 — 0.3 [panels (b) and 
(d)]. The largest real part of the complex eigenvalues A 
is shown by greyscale in the domain. It can be seen that 
increasing latency times reduce the domain of control. 
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FIG. 6: (Color online) Domain of control in the K-t plane for 
different latency times [panels (a) and (c): 5 = O.f , panels (b) 
and (d): 8 — 0.3]: The shaded areas indicate combinations 
of t and K, for which the largest real part of the complex 
eigenvalues A is negative and thus control is successful. The 
value of Re(A) is indicated by the greyscale (color online). 
The parameters of the unstable focus are chosen as u) = it in 
all panels and A = 0.5 in (a), (b) and A = 0.1 in (c), (d). 



For instance, the small range at a time delay of t = 1.5Tq 
in (a), where control is possible for S = 0.1, vanishes for 
6 = 0.3 in (b). Note that non-zero latency times lead 
to a loss of the symmetry [around r = (2n + 1)/2Tq for 
n = 1,2,...] of the domain of control (see also the case 
of zero latency as displayed in Fig. |3jl . 



V. LOW-PASS FILTERING 

It has been found that high frequency modulations of 
the control signal, due to additional high frequency com- 
ponents in the signal besides the main fre quen cy, can ren- 
der the TDAS control method unstable [33. As shown 
in that work, an additional low-pass filter included in the 
control loop can overcome this limitation, and UPOs can 
be stabilized [32L l33|. On the other hand, in electronic 
signal processing the finite response time of the circuit 
often imposes unavoidable low-pass filtering, and its ef- 
fect upon feedback control is not clear. In this section we 
will show that a low-pass filter changes the characteristic 
equation [see Eq. ©] of the fixed point, and shifts the 
minimal feedback gain to larger values. Note that low- 
pass filtering has been successfully used to stabilize USSs 
by generating a control force from the difference of the 
current state to its filtered counterpart [ill ]. 

The TDAS control force with an additional low-pass 
filter can be written as 



F(i) = -K 



x(t) — x(t — t) 
y{t) - V{t - r) 



(25) 




FIG. 7: Minimal feedback gain K vs. cut-off frequency a for 
A = 0.5 and u = n according to Eq. (1311) . The shaded area 
shows the domain of control. 



defined by 



x(t) = a / x{t')e~ a{t - t ' ) dt l 



(26) 



with the cut-off frequency a, and analogously for y(t). 
Equivalently, the convolution integrals can be replaced by 
two additional differential equations such that the origi- 
nal two-dimensional system becomes four-dimensional 

x(t) = \x{t) + uoy(t) - K[x(t) -x(t- r)] (27) 

y(t) = -uox{t)+Xy{t)-K[y{t)-y{t~T)] 

x(t) = —ax(t) + ax(t) 

y(t) = -ay(t) +ay(t). 

This system of differential equations yields a character- 
istic equation of the form 

± i(a + A)w = aK (l - e~ AT ) - (a + A)(A - A) (28) 

or, equivalently, using A — p + iq 

a(X-p) = aK [1 - e- pT cos(gr)] - q{q - uj) (29) 
a(u — q) = aKe~ pT sin^r) — Xq — pcu. 

Note that in the limit of large cut-off frequencies, i.e., 
a — > 00, Eqs. I)28|) and l)29|) reduce to the characteristic 
equations © and JSJ of Sec. [D] respectively. 

For further investigation of Eq. (|28|) we shall use the 
separation into real and imaginary parts (|29fl . Following 
the discussion of Sec. IIHI by considering even and odd 
multiples of tt as special values for qr, we obtain from 
the imaginary part of Eq. (|28fl in the case qr = (2n + l)ir 



2n+l 

T = 7T 



1 



(30) 



where x and y denote the filtered versions of x and y 



with n = 0, 1, 2, . . .. Inserting this into the real part of 
Eq. I|28|) gives an expression for the minimum value of 
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FIG. 8: (Color online) Largest real part of the complex eigen- 
values A vs. t for A = 0.5, u) = it, and K = 0.3 as given by 
Eq. 1281 . The solid, dashed, dotted, and dash-dotted curves 
correspond to a cut-off frequency of a = 1000, 10, 7 and 4, 
respectively. 



the feedback gain, for which the control method becomes 
successful for appropriately chosen r, 



K min (a) = — + 



A 



2 (a - A) 2 



(31) 



The dependence of the minimal feedback gain K m i n on 
the cut-off frequency a is depicted in Fig. [7| for A = 0.5 
and to = 7r. For large cut-off frequencies the minimum 
value tends to the result of Sec. IIIII [see Eq. (|llfl ]. Note 
that for finite a the minimal feedback gain is shifted to 
larger values compared to the case of the original TDAS 
control scheme. 

Figure [5] shows the largest real part of the eigenvalues 
A in dependence on the time delay r for fixed feedback 
gain K = 0.3 and various cut-off frequencies a = 1000, 
10, 7 and 4. For a large cut-off frequency, i.e., a — 1000, 
the curve is similar to the case without low-pass filter 
[see Fig. EJa)] indicating that the filter has only little ef- 
fect. For smaller a, however, filtering the control signal 
reduces the range of the time delay r, for which Re(A) 
becomes negative, eventually leading to a complete fail- 
ure of stabilization. Note that the notches are shifted to 
lower values of r for decreasing a. This effect can also 
be understood by Eq. (|30fl due to the additional factor 
(1 — A/a) < 1. The corresponding domains of control 
in the K-t plane are shifted to smaller values of r and 
distorted asymmetrically (Fig. |5J . 



puted the domain of stabilization of an unstable focus in 
the plane parametrized by feedback gain and time de- 
lay. Using the complex multivalued Lambert function, 
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FIG. 9: (Color online) Domain of control in the K-t plane 
for different cut-off frequencies [panels (a) and (c): a — 1000, 
panels (b) and (d): a — 7]: The shaded areas indicate com- 
binations of r and K , for which the largest real part of the 
complex eigenvalues A is negative and thus control is success- 
ful. The value of Re(A) is indicated by the greyscale (color 
online). The parameters of the unstable focus are chosen as 
uj — n in all panels and A = 0.5 in (a), (b) and A = 0.1 in (c), 
(d). 



we have derived analytically the main features of the 
stability domain by investigating the characteristic equa- 
tion of the fixed point. Below a minimum value of the 
feedback gain no control is possible. In the vicinity of 
this minimum value, the shape of the domain of control 
shows a square root dependence on the feedback gain. 
We find that no control is possible for time delays that 
are multiples of the uncontrolled eigenperiod of the sys- 
tem. Taking non-zero control loop latencies into account, 
we have shown that increasing latency times increase the 
minimum value of the feedback gain K m i n and reduce 
the domain of control substantially. Similarily, an addi- 
tional low-pass filter in the control loop causes a shift of 
K m in, as well. This suggests that filtering with a cut-off 
frequency a has a similar effect as a latency delay time 
or 1 . In fact, expanding 1)24(1 for small latency S in lowest 
order yields the same minimal feedback gain K'fnin as for 
low-pass filtering IpHJl. if Eq. (|3Tj|) is observed. 
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